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Abstract 

Thrichomys apereoides is an echimyid rodent which ranges in distribution from northeastern and cen¬ 
tral Brazil into Paraguay. Five subspecies are recognized, although each form is not well characterized 
and diagnosis is based primarily in pelage color variation. In this study we employed procedures from 
multivariate statistics to assess the systematic status of subspecies described from northeastern Brazil. 
The results of the craniometric analysis cannot be reconciled with the subspecies currently recognized 
for northeastern Brazil. Populations assigned to T. a. laiirentius and T. a. inermis form a continuum 
of variation in cranial size, although they differ in cranial shape from a population from the locality 
of Bodoco in the state of Pernambuco. The implications of these findings for the systematics of 
T. apereoides are discussed. 


Introduction 

Thrichomys apereoides (Lund, 1841) is an echimyid rodent commonly found in the savan¬ 
nas of eastern and central Brazil and Paraguay, often associated with rocky habitats 
(Mares et al. 1981; Alho 1982). This species also inhabits the dry caatingas of northeast¬ 
ern Brazil (Mares et al. 1981; Mares and Ojeda 1982), although it reaches higher densi¬ 
ties in mesic refuges. Thrichomys apereoides has been the subject of confusing 
nomenclatural and taxonomic assignements (Fetter 1973), and five subspecies are cur¬ 
rently recognized (Moojen 1952). The ranges of the subspecies are not well defined and 
the data available refer primarily to locality records. These are summarized in Moojen 
(1952), as follows: T. a. cunicularius known from the region of Rio Sao Francisco in Minas 
Gerais; T. a. pachyurus recorded from Paraguay and Mato Grosso and northern Sao Pau¬ 
lo; T. a. inermis described from Jacobina in the state of Bahia; T. a. apereoides ranging 
from western Minas Gerais to Goias; and T. a. laurentius occurring in the states of Ceara 
through Pernambuco. 

The subspecies of T. a. apereoides have been described primarily on the basis of pelage 
color variation (see Moojen 1952), and little is known about patterns of within and 
among-population variation in other character systems. In this study we employ multivari¬ 
ate morphometric procedures to examine variation in cranial metric traits in population 
samples of T. a. apereoides from northeastern Brazil which were allocated to the two sub¬ 
species, viz. T. a. laurentius and T. a. inermis, recognized for this area. The objectives of 
this study were (1) to assess the nature, extent, and magnitude of cranial differentiation in 
T. apereoides, and (2) to determine whether the pattern of population differentiation 
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based on multivariate analyses of cranial metric characters is consistent with the subspeci¬ 
fic structure established on the basis of pelage variation. 


Material and methods 

A total of 243 specimens of T apereoides available in the mammal collection of the Museu Nacional 
(Rio de Janeiro) was examined in this study. All specimens were classified to one of the seven age 
classes defined in Moojen et al. (1988) on the basis of tooth eruption and occlusal surface wear criteria. 
Only individuals belonging to age classes 5-7 were selected for the analysis of geographic variation be¬ 
cause they were considered adults by the criteria of Moojen et al. (1988). Statistically significant sexual 
size dimorphism has been detected for a few cranial traits in T. apereoides (Moojen et al. 1988), 
although in this study sexes were pooled to increase samples sizes. 

The specimens analyzed in this study represent samples collected at the following localities in north¬ 
eastern Brazil (Fig. 1): state of Ceara. Itapage (3°4r S, 39°35'W: n=13), Campos Sales (7°4'S, 
40°22'W: n = 7); state of Pernambuco. Bodoco (7°42'S, 39°53'W: n = ll), Triunfo (7°5r S, 38°8'W: 



A ITAPAGE 
B CAMPOS SALES 
BODOCO 
PRINCESAIZABEL 
TRIUNFO 
CARUARU 
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H SANTANA DO IPANEMA 
FEIRA DE SANTANA 
J PALMEIRAS 
* JACOBINA 


Fig. 1. Population samples of Thrichomys apereoides used in this study. Localities as indicated. The star 
denotes the type locality of Thrichomys apereoides inermis 





178 


A. C. Bandouk and S. F. Dos Reis 


n = 15), Camaru (8°14' S, 35°55' W: n = 9), Floresta (8°36' S, 38°44' W: n = 8); state of Parafba. Princesa 
Izabel (7°44' S, 37°59' W: n = 8); state of Alagoas. Santana do Ipanema (9°20' S, 37°16' W: n = 16); state 
of Bahia. Feira de Santana (12°15' S, 38°57' W: n = 9), Palmeiras (12°30' S, 41°34' W: n = 7). The popula¬ 
tion samples from the states of Ceara, Parafba, and Pernambuco were assigned to T. a. laurentius, and 
the two populations from Bahia were allocated to T. a. inermis (Moojen 1952). The population from 
Santana do Ipanema in the state of Alagoas does not fall into the ranges of either T a. laurentius or 
T. a. inermis and thus was not formally assigned to either of the subspecies. 

Thirteen quantitative cranial measurements defined in Smith and Patton (1988) in addition to two 
mandibular measurements were taken with electronic digital calipers, as follows: occipito-nasal length 
(ONL), basilar length (BAL), zygomatic breadth (ZB), mastoid breadth (MB), rostral length (RL), 
nasal length (NL), rostral width (RW), diastema (D), maxillary tooth row length (MTRL), palatal 
width (PW), bullar length (BUT), rostral depth (RD), cranial depth (CD), mandible length (ML), and 
mandible height (MH). 

Character variation in T. apereoides was analyzed by univariate and multivariate statistical proce¬ 
dures (Morrison 1976; Sokal and Rohlf 1981). Cranial size and shape variation among population 
samples of T apereoides was evaluated using principal components analysis (Morrison 1976), a multi¬ 
variate statistical procedure designed to summarize major patterns of covariation among metric traits. 
The values for all characters were transformed to logarithms to equalize variances and produce scale- 
invariant covariances that linearize allometric relationships (Bookstein et al. 1985). Principal compo¬ 
nents were extracted from the covariance matrix of log-transformed data. The principal component 
scores for the first three axes were plotted to assess patterns of cranial size and shape variation within 
and among populations of T. apereoides. 

All principal-component variate loadings are expressed as vector correlations (directional cosines; 
Morrison 1976), estimated for each character by its correlation with projection scores across indivi¬ 
duals. Vector plots graphically portray the principal directions of variation of each character within the 
corresponding principal component space. All statistical analysis were performed with SAS-PC (SAS 
Institute 1988). 


Results 

The statistics of variation in cranial dimensions in T apereoides are given in Table 1. Cra¬ 
nial dimensions do not show a clear pattern of variation in population samples from the 


Table 1. Mean {±SE) and analysis of variance for 15 cranial characters in populations of Thrichomys 
apereoides from ten localities in northeastern Brazil. 


Character 


Localities 



Itapage 

Campos Sales 

Bodoco 


X±SD 

X±SD 

X±SD 

Occipito-nasal length 

51.48 ±2.50 

49.1911.61 

49.7811.83 

Basilar length 

43.90 ± 2.28 

42.6611.20 

43.1611.76 

Zygomatic breadth 

25.71 ± 1.12 

24.6610.42 

25.0010.86 

Mastoid breadth 

19.8410.82 

19.1110.37 

19.4610.50 

Rostral length 

19.6411.22 

18.9010.95 

19.3511.11 

Nasal length 

17.9211.00 

17.05 10.87 

16.8910.91 

Rostral width 

9.5810.74 

8.9210.41 

10.5310.58 

Diastema 

9.9810.65 

9.33 10.44 

9.5710.48 

Maxillary tooth row length 

9.1510.25 

8.5210.19 

8.6110.42 

Palatal width 

2.9810.52 

3.1210.24 

3.43 10.30 

Bullar length 

11.2910.47 

11.4810.33 

11.6210.54 

Rostral depth 

10.6110.64 

9.9410.51 

11.6810.68 

Cranial depth 

13.2110.49 

13.1310.42 

14.6010.62 

Mandible length 

27.5711.43 

26.0610.71 

26.4110.86 

Mandible height 

11.8510.85 

11.3310.39 

11.5610.54 
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first four localities, whereas in the remaining populations there is an overall trend of de¬ 
creasing cranial dimensions from northern to southern localities. The decrease in mean 
values is more pronounced in the two southernmost populations from Feira de Santana 
and Palmeiras in the state of Bahia (Tab. 1). Univariate analysis of variation indicates sig¬ 
nificant heterogeneity among populations in all characters except nasal length (Tab. 1). 

Variation in craniometric measurements among populations of T apereoides was as¬ 
sessed by principal components analysis. The first principal component (PC) accounts for 


Table 1. Continued 


Character 


Localities 



Princesa Izabel 

Triunfo 

Caruaru 


X±SD 

X±SD 

X±SD 

Occipito-nasal length 

49.49 ±3.26 

52.3512.09 

51.3911.28 

Basilar length 

42.3013.45 

45.1012.14 

44.25 11.21 

Zygomatic breadth 

24.93 ±1.55 

25.9110.90 

26.0911.15 

Mastoid breadth 

19.5410.88 

20.2210.82 

20.1310.43 

Rostral length 

19.1411.66 

20.3711.07 

19.7510.60 

Nasal length 

17.0911.91 

18.2911.26 

17.5011.03 

Rostral width 

9.2410.66 

9.8210.79 

9.5710.69 

Diastema 

9.4010.85 

10.3310.45 

10.1610.39 

Maxillary tooth row length 

8.2810.75 

8.8010.21 

8.7910.31 

Palatal width 

2.9810.43 

3.3710.40 

3.4810.38 

Bullar length 

11.2710.47 

11.6310.64 

11.7510.37 

Rostral depth 

9.9910.78 

11.5610.59 

10.9210.65 

Cranial depth 

13.2710.78 

13.9910.63 

13.5310.52 

Mandible length 

26.3611.71 

27.7611.27 

27.6811.24 

Mandible height 

11.0610.80 

12.3010.66 

11.9010.64 


Table 1. Continued 


Character 


Localities 



Floresta 

Santana do Ipanema 

Feira de Santana 


X±SD 

X±SD 

X±SD 

Occipito-nasal length 

50.7712.02 

51.3812.53 

49.4712.09 

Basilar length 

43.85 11.66 

44.2812.40 

42.2411.77 

Zygomatic breadth 

25.4110.65 

25.2610.98 

24.5110.97 

Mastoid breadth 

20.0410.48 

19.7910.77 

18.8310.66 

Rostral length 

19.5910.83 

20.1111.25 

19.3811.26 

Nasal length 

17.6410.74 

17.7711.22 

17.4911.28 

Rostral width 

9.8910.89 

9.13 10.61 

9.6010.47 

Diastema 

10.0110.45 

10.15 10.92 

9.2010.55 

Maxillary tooth row length 

8.4110.41 

8.8010.30 

8.1410.30 

Palatal width 

3.5410.33 

3.2410.58 

3.2410.39 

Bullar length 

11.3910.56 

11.8710.57 

11.0810.53 

Rostral depth 

10.2910.29 

10.8410.67 

10.0310.62 

Cranial depth 

13.3810.46 

13.3210.58 

12.8810.60 

Mandible length 

27.7011.27 

27.2011.36 

25.5411.12 

Mandible height 

11.3810.67 

11.5510.95 

11.3910.63 
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Table 1. Continued 


Character 


Localities 


Palmeiras 


X±SD 

F 

P 

Occipito-nasal length 

47.35 ±1.65 

4.33 

0.0001 

Basilar length 

41.50 ±1.62 

2.93 

0.0042 

Zygomatic breadth 

23.66 ±0.74 

4.85 

0.0001 

Mastoid breadth 

18.26 ±0.69 

7.17 

0.0001 

Rostral length 

18.63 ±0.79 

2.26 

0.0243 

Nasal length 

16.75 ±0.53 

1.99 

0.0496 

Rostral width 

8.67 ±0.56 

6.09 

0.0001 

Diastema 

9.22 ±0.26 

5.08 

0.0001 

Maxillary tooth row length 

8.16 ±0.36 

8.45 

0.0001 

Palatal width 

3.24 ±0.32 

2.08 

0.0388 

Bullar length 

10.66 ± 0.41 

4.40 

0.0001 

Rostral depth 

9.60 ± 0.48 

12.96 

0.0001 

Cranial depth 

13.65 ±0.51 

7.89 

0.0001 

Mandible length 

25.03 ± 0.85 

5.72 

0.0001 

Mandible height 

10.03 ± 0.51 

6.40 

0.0001 


82.7% of the total variation, whereas PC-2 and PC-3 explain 4.3% and 3.0% of the varia¬ 
tion, respectively. Because all cranial characters were positively and significantly 
(P < 0.0001) correlated with the first pooled among-group principal component extracted 
from the covariance matrix of log-transformed character values, the first axis can be used 
as a general variable expressing overall cranial size. All populations overlap extensively 
in the space defined by the first two principal components (Fig. 2 A, B). On the other 
hand, the plot of PC-2 against PC-3 indicates the existence of two cranial shape morphol¬ 
ogies; one represented by the population of Bodoco (state of Pernambuco) and the other 
by a major cluster including all other populations (Fig. 3 A). The vector plot (Fig. 3 B) in¬ 
dicates that individuals in the population from Bodoco differ by having higher skulls and 
shorter rostra relative to individuals in the other populations. 

To further clarify the ordination of T. apereoides we performed canonical discriminant 
function analysis (Morrison 1976). The pattern of ordination is similar to that uncovered 
with principal components analysis (Fig. 4 A), in that the population from Bodoco is fully 
discriminated. Ordination of the remaining populations is more informative in the canoni¬ 
cal discriminant analysis. The populations from Itapage, Campos Sales, Triunfo, Caruaru, 
Floresta, Princesa Izabel, Santana do Ipanema, Feira de Santana, and Palmeiras are dis¬ 
tributed along an axis (CV-2) of increasing cranial size starting with the population from 
Itapage and ending with the population from Palmeiras (Fig. 4 A). Canonical variate 2 is 
clearly an axis of size since all cranial characters have positive correlations with this vec¬ 
tor (Fig. 4B). 


Discussion 

The statistical analysis of cranial metric traits in T. apereoides in northeastern Brazil re¬ 
vealed the existence of two groups of populations which lie in discrete regions in multi¬ 
variate character space and differ in cranial shape. One major group is widespread in 
distribution and includes population samples assigned to T. a. laurentius and T. a. inermis 
from the localities of Itapage, Campos Sales, Triunfo, Caruaru, Floresta, Princesa Iza- 
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Correlation with PCI 


Fig. 2. Principal components analysis for ten populations of Thrichomys apereoides. (A) Bivariate plots 
of projection scores on principal components 1 and 2. Letters at the tip of polygons indicate sample lo¬ 
cality and letters at the middle of polygons indicate sample centroids. (B) Vectors portraying the princi¬ 
pal directions of character variation in the plane of the first two principal components 
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Correlation with PC2 


Fig. 3. Principal components analysis for ten populations of Thrichomys apereoides. (A) Bivariate plots 
of projection scores on principal components 2 and 3. Letters at the tip of polygons indicate sample lo¬ 
cality and letters in the middle of polygons indicate sample centroids. (B) Vectors portraying the princi¬ 
pal directions of character variation in the plane of the principal components 2 and 3 
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Correlation with CVl 


Fig. 4. Canonical discriminant analysis for ten populations of Thrichomys apereoides. (A) Bivariate 
plots of projection scores on canonical variates 1 and 2. Letters at the tip of polygons indicate sample 
locality and letters at the middle of polygons indicate sample centroids. (B) Vectors portraying the prin¬ 
cipal directions of character variation in the plane of the first two canonical variates 
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bel, Santana do Ipanema, Feira de Santana, and Palmeiras. Variation in this group is 
structured in a dine of decreasing cranial dimensions from northern to southern local¬ 
ities. The second group is represented by a single population from Bodoco, which was 
allocated to T. a. laurentius, and does not fall on the dine of decreasing cranial dimen¬ 
sions. 

The morphologic evidence presented here for the patterns of variation in 
T. apereoides can be interpreted according to current views on the structure of subspeci¬ 
fic variation. Barrowclough (1982) and Thorpe (1987) have stressed that evidence for 
the recognition of subspecies should not be based on continuous clinal variation, but 
that populations should occupy discrete regions in character space. The nature of cranial 
differences among populations of T. apereoides is also relevant to the definition of infra¬ 
specific differentiation. The two components of form, i. e., size and shape have been at¬ 
tached different weights in taxonomy and systematics, with the latter component (shape) 
being held as a better predictor of genetic variation and relatedness and phytogeny 
(Gould and Johnston 1972; Thorpe 1983; Rohlf and Bookstein 1987). Recently, Pat¬ 
ton and Brylsky (1987) demonstrated that size is phenotypically plastic in populations 
of the pocket gopher Thomomys bottae, and shifts in cranial size are a response to en¬ 
vironmental (nutritional) quality. On the other hand, cranial shape variation was found 
to be correlated with genic differences in recognized subspecies of T. bottae. Smith and 
Patton (1988) elaborated further on the relation between cranial size and shape and 
ecologic and historical (phylogenetic) components of differentiation, and suggested that 
independent evolutionary units (= subspecies, sensu Smith and Patton, 1988) should 
be recognized for geographic units showing concordant patterns of cranial shape and 
genetic variation. 

Based on current view of the structure of subspecific variation (cf. Barrowclough 
1982; Patton and Brylski 1987; Thorpe 1987; Patton and Smith, 1989, 1990), the cranio- 
metric data suggest the existence of two geographic units based on cranial shape: one re¬ 
presented by the population from Bodoco and the other including populations from 
Itapage, Campos Sales, Triunfo, Caruaru, Floresta, Princesa Izabel, Santana do Ipanema, 
Feira de Santana, and Palmeiras. This result is not consistent with Moojen’s (1952) sub¬ 
specific arrangement of T. apereoides in northeastern Brazil and raises two points. First, 
T. a. laurentius, which according to Moojen (1952) ranges in distribution from the state of 
Ceara to the state of Pernambuco, may not be monophyletic but rather be composed of 
two geographic units. Second, the populations from Feira de Santana and Palmeiras in the 
state of Bahia which were allocated to T. a. inermis, were shown in fact to be continuously 
connected with other populations assigned to T. a. laurentius. Our samples from the state 
of Bahia are not from Jacobina, the type locality of T. a. inermis, and thus definition of 
the status of this subspecies can only be established when material from this locality be¬ 
comes available. Additional sampling will also be necessary to determine the geographic 
extent of the cranial morphology represented by the population from Bodoco. A finer 
geographic sampling of the cranial morphology coupled with surveys of genic and chro¬ 
mosomal variation will be needed to better understand the structure of variation and evo¬ 
lution in T. apereoides. 
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Zusammenfassung 

Kraniometrische Variation and Unterartgliedening bei Thrichomys apereoides in Nordost-Brasilien 
(Rodentia: Echimyidae) 

Das Verbreitungsgebiet der Stachelratte Thrichomys apereoides reicht von Nordost- iiber Zentral-Bra- 
silien bis nach Paraguay. Funf Unterarten werden gegenwartig anerkannt, auch wenn diese Formen 
nicht deutlich charakterisiert sind und ihre Diagnosen meist auf Farbungsmerkmalen basieren. Fur die 
vorliegende Studie wurden Schadel von 243 Exemplaren vermessen und die Werte einer multivariaten 
statistischen Analyse unterzogen. Die Resultate stimmen nicht mit der gegenwartigen Unterartglieder- 
ung uberein. Populationen, die traditionell als T. a. laurentius und T. a. inermis bestimmt wurden, bilden 
in bezug auf ihre SchadelmaBe ein Kontinuum, unterscheiden sich andererseits aber in der Schadelform 
von einer Population aus Bodoco im Bundesstaat Pernambuco. Innerhalb T. apereoides gibt es offenbar 
zwei gut unterscheidbare Populationen, deren systematischer Status diskutiert wird. 
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